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ABSTRACT 

We describe a new method for measuring galaxy magnification due to weak gravitational lensing. 
Our method makes use of a tight scaling relation between galaxy properties that are modified by grav- 
itational lensing, such as apparent size, and other properties that are not, such as surface brightness. 
In particular, we use a version of the well-known fundamental plane relation for early type galax- 
ies. This modified "photometric fundamental plane" replaces velocity dispersions with photometric 
galaxy properties, thus obviating the need for spectroscopic data. We present the first detection of 
magnification using this method by applying it to photometric catalogs from the Sloan Digital Sky 
Survey. This analysis shows that the derived magnification signal is comparable to that available 
from conventional methods using gravitational shear. We suppress the dominant sources of system- 
atic error and discuss modest improvements that may allow this method to equal or even surpass the 
signal-to-noise achievable with shear. Moreover, some of the dominant sources of systematic error 
are substantially different from those of shear-based techniques. Thus, combining the two techniques 
addresses the major weaknesses of each and provides a substantial improvement over either method 
used in isolation. With this new technique, magnification becomes a necessary measurement tool for 
the coming era of large ground-based surveys intending to measure gravitational lensing. 

Subject headings: cosmology: observations — gravitational lensing: weak — methods: observational 



1. INTRODUCTION 

The modern cosmological concordance model has been 
spectacularly successful. This success has come at a 
price, however: cosmologists must postulate a universe 
in which the vast majority of the content, in the form 
of dark matter and dark energy, is inaccessible to direct 
observation. 

Gravitational lensing, and weak lensing in particular, 
provide our best window onto the dark universe. Because 
of this, the astronomical community is investing heavily 
in current and future imaging surveys, both ground- and 
space-based, designed at least in part around weak lens- 
ing science, e.g., the Dark Energy Survey (DES), the 
Panoramic Survey Telescope and Rapid Response Sys- 
tem (Pan-STARRS), the Hyper Suprime Cam for the 
Subaru telescope, the Large Synoptic Survey Telescope 
(LSST), and Euclid. 

Lensing measurements have already played a signifi- 
cant role in astrophysics in the last two decades over a 
range of scales and physical regimes. Weak lensing mea- 



surements have characterized the agg 
galaxies' dark matter haloes, (e.g., 
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120091 : iLeauthaud et al.1 12010ft and, with recent co smic 
shear detections (jFu et alJl2008t iMassey et al.ll2007f ). di- 
rectly measured the clustering of matter on cosmological 
scales. 

These measurements are currently made almost ex- 
clusively by studying spatially-correlated distortions in 
the ellipticities of background galaxies due to the shear 
component of the gravitational lensing distortion. In 
the weak-lensing regime, the ellipticities induced by lens- 
ing are small (~ 1%) compa red to the range of intrin sic 
galaxy ellipticities (~ 30%) (IMandelbaum et al.|[2005T) . 

The lensing distortion also magnifies background 
sources, but the intrinsic variance in the distribu- 
tion of galaxy sizes and luminosities — those proper- 
ties perturbed by magnification — is much larger than 
that of galaxy shapes. Magnification measurements to 
date have necessarily had to average over much larger 
galaxy samples to obtain sign al-to- noise (S/N) equiv- 
alent to shear measurements (|Hildebrandt et all 120111; 
Bauer et"aTl 120111 : [Menard et al.l l2010t iScranton et afl 
2005). 

In this Letter, we make use of 
scaling relation — a photometry-only 
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the dark matter profiles of large ga laxy clusters , both 
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well-known fundamental plane for early type galaxies 
(jDiorgovski fc Davisiri987l : lDressler et al.lll987l )— to sub- 
stantially narrow the intrinsic distribution of sizes for a 
set of background source galaxies. This makes it possible 
to measure the weak lensing magnification signal around 
a set of foreground galaxy lense s with higher S/N tha n 
was previously thought possible (Rozo & Schmidt 2010). 
A thorough description of the methodology and imple- 
mentation of this analysis will be presented in a forth- 
coming paper (hereafter Paper II) . Here, we describe the 
core concepts, present a detailed outline of the method, 
and demonstrate its effectiveness. 
In section [2 we describe the scaling relation used in 
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this analysis. In section [31 we present a proof-of-concept 
measurement of weak lensing magnification using this 
technique and control for the most important systematic 
biases. We conclude in section [4] with a brief discussion 
of the ways this method might be improved upon, with 
an eye toward extracting a lensing signal from magnifi- 
cation that equals or even exceeds the S/N obtainable 
from shear-based techniques. Throughout, we assume a 
ACDM cosmology with M = 0.274, D, A = 0.726, and 
H = 100 km s" 1 Mpc" 1 . 

2. THE PHOTOMETRIC FUNDAMENTAL PLANE 

The fundamental plane (FP) is in many ways an ideal 
tool for measuring magnification. It is an observed corre- 
lation between galaxy effective radius (Re), which is mag- 
nified by gravitational lensing, and two galaxy properties 
which are unaltered by lensing: galaxy surface brightness 
(/u) and the stellar velocity disper sion (a). The intrinsic 
scatter in the FP is ~ 0.08 dex (| Jorgensen et al.l fl996l : 
iBernardi et all [20031 ) . or 20%. Thus the FP makes it 
possible to predict the intrinsic value of R e from obser- 
vations of [i and er, which can then be compared with the 
observed values of R e to measure magnification. 

The FP was in fact propose d as a tool for this purpose 
bv lBertin fc Lomba rdi (2006), but to our knowledge has 
never been used as such due to a critical flaw. Placing 
galaxies on the FP requires a measurements. Even with 
the tight scatter in the FP, a statistically viable measure- 
ment would require high-resolution spectroscopic mea- 
surements for millions of galaxies. 

Identifying a purely photometric analog to the FP with 
comparable scatter would solve this p roblem. Such a 
relation has already been identified by Graham (2002), 
where the concentration of the galaxy light profile fills the 
role normally served by a. This works in part because 
concentration and velocity dispersion are both strongly 
correlated with galaxy mass, and in part because at fixed 
mass, galaxies with more concentrated mass profiles have 
higher velocity dispersions. We will explore the rela- 
tion between the spectroscopic fundamental plane and 
the photometric relation deployed here more fully in Pa- 
per II. 

2.1. Background Sources 

To define a photoFP for this work, we use a sample 
of galaxies drawn from the Sloan Digital Sky Survey 
III (S DSS-III) Eighth Data Release (DR8, lAihara et all 
1201 If ). We limit the sample to resolved sources that meet 
basic quality cuts (e.g., are not saturated). For these, 
we estimate photometric redshifts (photo-z's) based on 
the S PSS uqriz photometry using the public code ZE- 
BRA (jFeldmann et al.ll2006f l run with the default tem- 
plates, allowing interpolation between the standard tem- 
plates without template optimization. To select a sam- 
ple of early type background sources that should lie on 
the photoFP, we exclude the ~ 2/3 of the galaxies with 
best-fitting templates inconsistent with that of a passive 
stellar population. The sample selection for background 
sources will be described in greater detail in Paper II. 

The SDSS photometric pipeline does not measure 
Sersic index. Here, we substitute for n the SDSS pet- 
rosian concentration C = R90/R50, defined as the ratio 
of the radii containing 90% and 50% of the Petrosian flux 




-10 12 
f(* C) 

Figure 1. The photometric fundamental plane for our source sam- 
ple of 8.4 million galaxies, shown edge-on. log R e is fit as a function 
of effective surface brightness (fi) and concentration (log C) sepa- 
rately in redshift bins of width Az = 0.01. Gray points show a 
random subset of 100,000 galaxies from the source catalog, while 
the solid line shows the one-to-one relation. Contours enclose the 
0.5(T, la, 1.5it, 2a, 2.5(7, and 3a boundaries of the 2D distribution 
for the full source catalog. The inset shows the distribution of resid- 
uals in logi? e from the photoFP fits, which has width a = 0.153 
dex. 

(e.g.. lShimasaku et al.|[200l . All reported quantities are 
measured in the r band. 
We fit a photoFP of the form 

logi? e = a/i + /31ogC + 7, (1) 

where R e is the half-light radius of the best-fit de Vau- 
couleurs light profile converted into physical units us- 
ing the ZEBRA photo-z, /i is the mean de Vaucouleurs 
surface brightness within R e , and a, (3, and 7 are free 
parameters. To avoid errors resulting from a redshift- 
dependent selection function, evolution in the photoFP, 
and if-corrections to the radii due to the fact that 
the morphological measurements are all made in the 
observed-frame r band, we divide our galaxy sample into 
redshift bins with width Az = 0.01 and fit the photoFP 
separately in each bin. The best-fit coefficients are cho- 
sen to minimize the dispersion in effective radius at fixed 
/i and logC, taking into account only the errors in R e . 

Figure Q] shows an edge-on view of the photoFP for 
our source sample. The dispersion around the photoFP 
in the direction of effective radius is 0.15 dex, or 35%. 

2.2. Magnification using the photoFP 

A line-of-sight matter overdensity at lens redshift zi 
will produce an image convergence k of amplitude: 



■^crit 

where £ is the projected surface density on the sky at 
zi and S C rit is the characteristic surface density of mat- 
ter required for lensing. S cr ;t is defined by the lensing 
geometry, such that 

c 2 d s 
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Figure 2. The projected correlation function for sources around 
lenses with 0.20 < z; < 0.25. Colors indicate different bins in z s . 
Solid lines show fits to the data. At small separations, a large 
fraction of the galaxies in nearby z bins are likely scattered in from 
Zi through photo-z errors. See text for details. 



where di, d s , and di s are the angular diameter distances 
from the observer to the lens, from the observer to the 
source, and from the lens to the source, respectively. The 
factor of (1 + z{) 2 arises from our use of comoving coor- 
dinates. 

The lensing convergence re-scales the light profile, in 
the limit of very weak lensing, by a factor of (1 + k). The 
radius and luminosity increase, but as the light profile is 
simply rescaled, the concentration is left unchanged. In 
the presence of the scaling relation described above, this 
implies an estimator k of: 



log (1 + ft) = Alogi? e 



logi? e - (o^t + /3 log C + 7). 



(4) 



If the errors in the observables are uncorrelated, the 
variance in our estimator k is just the variance in the 
photoFP in the direction of R e . We extract a galaxy- 
galaxy lensing signal by cross-correlating this estimator 
with a population of foreground lenses. 

3. A MAGNIFICATION MEASUREMENT 

3.1. Lens Sample 

The le ns sample is selected f rom the NYU Value-Added 
Catalog (|Blanton et all 120051) version o f the SPSS Data 
Relea se 7 (DR7) spectroscopic survey ()Abazaiian et al.l 
2009), using only Luminous Red Galaxy Sample tar- 
gets (LRGs, lEisenstein et al.l I2003D . In order to com- 
pare with the results of iMandelbaum et al.1 (j2008al ). we 
limit the sample to massive galaxies with absolute r- 
band magnitudes —21.5 > Mo.o r > —22.6 and redshifts 
0.15 < z < 0.35. The magnitudes are fc-corrected and 
evoluti on corrected to z = 0.0 as in Mandcl baum et al.l 
(|2006al hereafter M+06). Finally, to exclude satellite 
galaxies that are not at the centers of their dark matter 
haloes, we remove galaxies with brighter nearby LRGs, 
again following M+06 . This gives a sample of ~ 55, 000 
lenses that have comparable properties to the combined 
LRG sample of M+06. 



3.2. Correcting Biases due to Photometric Redshift 
Errors 

In the presence of photo-z errors, the overdensity of 
sources clustered near a lens will produce an excess of 
galaxies with incorrect photo-z (z p ) along the line of sight 
to the lens. As a result, when we average Alogi? e over 
the foreground or background source galaxies, we sys- 
tematically mis-estimate the residuals from the plane as- 
sociated with a lens due to the 'shadow' cast by photo-z 
errors. 

To deal with this bias, we will calculate the magnitude 
of this spurious signal directly from the data, and sub- 
tract it from our measured signal. We must first estimate 
the error in A log R e induced by a galaxy being assigned 
the wrong z p (A log i?g 1T ), then calculate what fraction 
/; of the galaxies at each z p have been scattered in from 
z\. In these terms, the observed mean photoFP residual 
is: 



A log R° = (1 — fi) log (! + «) + /,Alogi?° 



(5) 



where k is the true convergence. 

A log R°" can be estimated by assuming that the 
galaxy lies on the photoFP at Zi but is incorrectly as- 
signed to z p . The inferred effective radius of a galaxy 
with true redshift zi that is mistakenly assigned to z p will 
be off by a factor of d s (z p ) /d s (zi). The surface bright- 
ness dimming correction will be similarly incorrect, with 
Hp = w - 10 log [(1 + z p ) I (1 + z{)]. Finally, the photoFP 
fits differ between redshift bins. A galaxy with an in- 
correct photo-z will therefore lie off the photoFP in its 
assigned redshift bin by 



A log = log 



d a (z p )RP(ii p ,C) 
d s (zi) R l e (n h C) 



(6) 



The expressions RP e (fi p ,C) and R l e (ni,C) are the radii 
that would be predicted by the photoFP for that galaxy's 
surface brightness and concentration in the bins corre- 
sponding to z p and zi, respectively. 

The quantity /; can be estimated by cross-correlating 
the positions of sources at z p with lenses at zi. We as- 
sume that the positions of galaxies in widely separated 
redshift bins are uncorrelated and that any observed ex- 
cess of sources far behind a lens is due to scattering from 
zi. This means that 



fl = 



wg(9) 

\ + w a (ey 



(7) 



where wu (9) is the angular cross-correlation between the 
positions of sources at Zi and lenses at z\. A cross- 
correlation signal of this form can also be produced 
by the booste d number counts o f magnified background 
sources (e.g., Uain h Limall2011r ) but that effect is too 
weak to detect with a lens sample of this size. 

The cross-correlations for 0.20 < z\ < 0.25 with a range 
of z s bins are shown in figure [5] We fit an angular cor- 
relation function of the form 



1 + Wtt{0) 



An 

00.8 



■B* 



(8) 



where An and Bn are free parameters. The choice of 
power law index is motiva ted by the angular correlation 
function measurements of I Wake et al.l (|2011l ). which are 
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in agreement with our observed wu ■ Incorrectly estimat- 
ing the true mean density of galaxies at z p will cause Bu 
to deviate from unity, as is observed. We remove the 
effects of this uncertainty when calculating /; by setting 
Bu = 1. Sources with /; > 0.20 (above the black hor- 
izontal line in Figure 2) are excluded from the lensing 
measurement, while sources with /; < 0.20 are corrected 
using equation [5] 

In addition to the effects of galaxy clustering on pho- 
tometric redshift errors, a mean offset between the true 
and photometric redshifts in a z p bin will cause an in- 
correct estimation of the critical density £ C rit for all of 
the galaxies in that bin. This error depends on the dis- 
tri bution of foregroun d lens re dshifts. Using the method 
of Mandclbau m et all (|2008bf) , we estimate the effect of 
a mean shift in our photo-z's on the signal of no more 
than 10%. This uncertainty is small relative to the other 
corrections discussed here, so we defer this calculation to 
Paper II. 

3.3. Sky Proximity Bias Correction 

The SDSS photometric pipeline produces known sky 
subtraction proximity effects, where the photometry of 
objects ne ar bright stars or ga laxies is systematically bi- 
ased (c.f. lAihara et all 120111 ). This may induce a sys- 
tematic bias in the estimated radii, surface brightnesses, 
and concentrations that contaminates the lensing signal. 
Sky subtraction effects cannot distinguish between fore- 
ground and background galaxies (with respect to the 
bright lens), so this proximity bias can in principle be 
estimated from the photoFP residuals for galaxies in the 
foreground of the lenses, which are unaffected by lensing. 

Figure[3]shows the average deviation from the photoFP 
as a function of source-lens angular separation for both 
foreground and background sources. The sky proximity 
bias systematically induces a reduction in effective radius 
relative to the photoFP trend. The lensing signal is thus 
the difference between the background and foreground 
photoFP deviations at each angular separation. Of note 
is the fact that our empirical sky correction extends be- 
yond the size of the SDSS sky subtraction box, which 
is ~ 100"; this is a result of galaxy clustering. Each of 
the bright objects used as a lens will tend to be associ- 
ated with a galaxy overdensity on the sky. This excess 
will also impact the sky correction, even in neighboring 
sky subtraction cells, so the angular scale of the resulting 
correction will be set by the galaxy correlation function. 

The low redshift of the lens sample and the poor qual- 
ity of the photo-z's (which preferentially scatter higher- z 
sources to lower z) result in a large fraction of source 
galaxies near lenses with photometric redshifts z p < z\ 
that are actually at z > zi . This means that a foreground 
sample of sources with photo-z's will be contaminated by 
objects from higher z. The cut on described above re- 
moves many such contaminating galaxies, at the cost of 
dramatically reducing the signal-to-noise ratio of the sky 
proximity bias estimate. This remains the major source 
of uncertainty in this measurement. 

As a check against this effect, we also show the devia- 
tion from the photoFP trend of those foreground sources 
with spectroscopic redshifts. Any large bias to this sky 
subtraction estimate resulting from imperfect photo-z's 
should produce a substantial difference between the spec- 
troscopic and photometric foreground estimates; this is 
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Figure 3. The raw magnification signal around the galaxy lenses 
(filled black circles) compared with the sky proximity bias mea- 
sured from foreground sources. The red crosses show the estimated 
sky subtraction effect using sources with photo-z's; the blue trian- 
gles show the same estimate, but using those foreground galaxies 
with spectroscopic redshifts. 

not observed. 

3.4. Halo Mass Profile 

After controlling for the systematic errors described 
above, we calculate the line-of-sight surface matter den- 
sity S by weighting each lens-background source pair by 
the critical density for lensing, £ cr it (z s , zi). We bin this 
density by physical separation in the lens plane. Our 
results are shown in figure SI along with existing mea- 
surements from M+06 for a similar lens population. 

4. DISCUSSION: THE WAY FORWARD 

The magnification signal demonstrated above, while 
many times stronger than previous magnification mea- 
surements, is still somewhat noisier than the shear signal 
for a comparable sample. This is because the convergence 
dispersion resulting from the measured photoFP width is 
35% (1.8 times larger than the intrinsic shear dispersion 
of 20%) and because we have only used the third of the 
source sample consistent with early-type SEDs. 

If the fundamental achievable limit for this technique 
is the intrinsic scatter in the spectroscopic fundamental 
plane, then the average magnification S/N for an early- 
type galaxy is the same as in shear; a comparable photo- 
metric Tully-Fisher relation for late-type galaxies would 
bring us to the point where magnification and shear pro- 
vide comparable information. And any improvement in 
our understanding of galaxy evolution and dynamics that 
further diminishes the scatter in these scaling relations 
will boost the magnification signal beyond that available 
for shear measurement. 

Perhaps just as valuable, magnification by this method 
is not sensitive to the same systematic biases that chal- 
lenge upcoming shear measurements. For instance, the 
intrinsic galaxy alignment signal on large scales should 
not affect galaxy sizes, concentrations, and mean sur- 
face brightnesses in same the manner in which it affects 
shapes. We expect that this technique will also prove 
useful in extracting and removing instrumental sytemat- 
ics, such as those arising from variations in the telescope 
point-spread function, and will investigate this prospect 
in a subsequent paper. 
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Figure 4. Red points: S from this work. Black points: AS from M+06 measured using shear. That measurement used a smaller lens 
sample than we consider here, so we have reduced those error bars to allow for a fair comparison with this sample. The solid black line is 
the best-fit AS profile from M+06. The shaded red region shows the corresponding S profile (with 68% confidence interval) derived from 
the M+06 data. 
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